February 1, 2008 



LBNL-42429 
UCB-PTH-98/53 
hep-ph/9810468| 



Tau Appearance In 
Atmospheric Neutrino Interactions!] 

Lawrence J. Hall and Hitoshi Murayama 

Department of Physics 
University of California, Berkeley, California 94720 

Theoretical Physics Group 
Ernest Orlando Lawrence Berkeley National Laboratory 
University of California, Berkeley, California 94720 

Abstract 

If the correct interpretation of the Super-Kamiokande atmospheric 
neutrino data is — > v T oscillation, the contained data sample should 
already have more than 10 r appearance events. We study the chal- 
lenging task of detecting the r, focussing on the decay chain — ► 

— ► vr^vr in events with quasi-elastic r production. The background 
level, which is currently quite uncertain because of a lack of relevant 
neutral current data, can be measured by the near detector in the K2K 
experiment. Our estimates of the background suggest that it may be 
possible to detect r appearance in Super-Kamiokande with 5-10 years 
of running. 
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1 Introduction 



Recent data from the Super-Kamiokande (SK) detector has provided 
confirmation of the atmospheric neutrino anomaly, previously reported by 
Kamiokande @ and 1MB f| collaborations. The Soudan-II f| and MACRO 
H Collaborations have also reported anomalies which support the Super- 
Kamiokande observations. Apparently about 40% of the produced in cos- 
mic ray showers in the atmosphere do not reach the underground detectors. 
While calculations of the flux contains uncertainties, these calculations 
can be checked by comparing with the observed rate of muon and u e induced 
events. Neutrino oscillations provide a plausible depletion mechanism which 
has been strengthened by recent Super-Kamiokande data showing that the 
amount of depletion depends on the azimuthal angle. 

In this letter, we argue that, if the oscillation is to u T , then the Super- 
Kamiokande detector could detect charged current r appearance events via 
the chain u T — > r~ — > p~ — > 7r~7r°. The signal event rate is small because only 
a small fraction of u T have sufficient energy to produce r. Furthermore, there 
is considerable uncertainty in the background from p production via other 
neutrino induced processes. However, we estimate the signal to background 
ratio to be above unity, and the signal has a strong up-down asymmetry. A 
multi-year search at Super-Kamiokande with an efficiency of 50% could lead 
to a 5cr signal. 

The atmospheric could oscillate into u e , v T or u s , an SU(2) singlet neu- 
tral fermion. Fits to the atmospheric data, in models with three generations 
and a single Am 2 relevant to the atmospheric oscillation, show that at most 
about a third of the oscillations can be to v e , so that at least two thirds are 
to v T 0. For values of Am 2 above 2 x 10~ 3 eV 2 , the CHOOZ experiment 
constrains the v e fraction to be even smaller J7]. Within the minimal sce- 
nario, the oscillations are dominantly to v T . While oscillations to v s remain 
a possibility, it is theoretically less favoredJJ] 

How can the preferred — > v T interpretation be verified? A long baseline 
experiment with low neutrino beam energy (e.g., KEK to SK) would see 

tit requires the addition of an exotic fermion, neutral under the known gauge inter- 
actions. Many theories contain such fermions; however, because of their neutrality, they 
typically get very large masses, thereby inducing small masses for v e ,n,r by the see-saw 
mechanism. On the other hand, for to oscillate into v s , the mass of v s must certainly 
be less than 250 keV. Understanding such a light mass for a singlet state is possible, but 
requires more elaborate theoretical structures. 
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disappearance but not r appearance. For higher energy neutrino beams (e.g., 
MINOS and ICARUS) the more direct signal of r appearance is possible. 
In this letter we argue that a r appearance signal, with a large up-down 
asymmetry, may be possible at Super-Kamiokande. We caution that the 
event rates are low, and that a significant p detection efficiency is necessary. 
It may be possible to see this signal even if Am 2 < 10~ 3 eV 2 , where the 
planned accelerator based experiments become more difficult. 



2 t Appearance 

In this section, we give numerical results for r-appearance cross sections, 
from both quasi-elastic (QE) processes 

v T n -> r~p, (1) 
u T p — > r + n, (2) 

and deep inelastic scattering (DIS) processes 

u T p -> t~X, (3) 

v r n —> r~X, (4) 

v T V -> r + X, (5) 

v T n -> t + X. (6) 

We have not calculated the additional contributions from resonance produc- 
tion. Therefore, the sum of the QE and DIS contributions gives a conservative 
estimate of the total r production cross section. 

The QE cross sectionsare evaluated using the formula in Ref. [§], with the 
second-class currents dropped, PCAC and dipole form factors assumed, and 
taking rriA — 1-0 GeV as suggested by experiments ||. The results are shown 
in Fig. [I] together with the corresponding quasi-elastic cross sections for muon 
(anti-) neutrino induced events. We quote cross sections for r production at 
E v = 5 GeV: 

°QE(v T n -> r~p) = 3.5 fb, (7) 
a QE (v T p^T + n) = 2.0 fb. (8) 

The DIS cross sections are evaluated using the MRS-A parton luminosi- 
ties, specifically extended to the low-Q 2 regime [10]. The parton-level cross 
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Figure 1: The cross sections of r production from v T charged current reac- 
tions, for quasi-elastic (QE) events as described in the text. The \i production 
cross sections are shown for comparison. 



sections, however, are the leading-order calculations fLl |. Cross sections for 
all four processes, (^]||,|5],|6|) are shown in Fig. ||.| 

The QE cross sections are combined with calculations of the atmospheric 
neutrino flux [12 for a 414 day exposure at SuperKamiokande (21.1 kton), 
assuming no nuclear effects. The total number of electron events calculated 
is 106.7 for the multi-GeV sample, p e > 1.33 GeV. According to the MC 
calculations based on the same atmospheric neutrino flux [12] shown in the 
2nd paper in [jl|, the Super-Kamiokande collaboration expect 182.7 e-like 
events, with 38.3% from the QE process. Therefore, we overestimate the QE 
event rate, due to our lack of knowledge of the detection efficiency and the 
nuclear shadowing effect, and we should multiply our calculated results by a 
reduction factor of 0.66. 

To confirm this reduction factor, we further compared the multi-GeV fi- 

all of our discussions, we took the r threshold from nucleons at rest. In reality, 
however, the Fermi motion of nucleons inside oxygen nuclei lowers the threshold and hence 
the expected signal rate is higher. We thank G. Shapiro on this point. A precise estimate 
of this effect, however, requires a nuclear model and is beyond the scope of this paper. It 
is worth emphasizing though that we estimate the signal rate conservatively only using 
QE process with fixed production threshold. 
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Figure 2: The cross sections of r production from v T charged current reactions 
for deep-inelastic (DIS) events, normalized by G 2 F m^E v j^ . The rising solid 
curves [y T CC) are for v T n, u T p, u T p, and v T n from below. The neutral current 
[y NC) processes are given in dashed lines, for z/p(n), up, vn from below; the 
first two are accidentally close and are not distinguishable in the plot. The 
CC production cross sections are shown for comparison, for u^n, v^ip^p, 
and Vyji from above, where the middle two are accidentally close and are not 
distinguishable in the plot. 
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like events to our calculations. This is slightly more complicated because 
there are both fully contained (FC) and partially contained (PC) events. 
We find that there should be 250.0 QE events. Their MC shows 229.0 FC 
events, with a 54.1% QE fraction, and 287.7 PC events with 17.9% QE: 
in total, 175.5 QE events. The same reduction factor of 0.66 explains the 
descrepancy (actually about 0.7), and hence we adopt this factor to scale 
between our calculated rates and those from the Super-Kamiokande Monte 
Carlo. 

The Super-Kamiokande Collaboration observed 218 e-like multi-GeV events 
in 414 days. This is larger than the 182.7 events expected from their Monte 
Carlo calculation by a factor of 1.19. We assume that this has nothing to 
do with oscillations, but is due to the uncertainties in the flux calculation. 
We therefore assume that all unoscillated fluxes, produced in the cosmic ray 
showers, are a factor 1.19 larger than calculated in []12[ , which increases our 
reduction factor from 0.66 to 0.77. 

We studied the momentum distribution of e and \i events after rescaling 
the normalization as discussed, for a 414 day exposure at SuperKamiokande 
(21.1 kton). This distribution for [i events is shown in Fig. |3|(a). Compared 
to the reported data, the fraction of QE events in the one-ring event samples 
becomes smaller for higher momentum as expected on physical grounds. The 
other events correspond to the production of a few additional soft particles 
below the detection threshold. We will come back to this issue later. 

Including v il -v T oscillations for the "test point" used by Super-Kamiokande 



Am 2 = 0.005 eV 2 , sin 2 9 = 1, (9) 

we obtain the dashed histogram in Fig. |3|(a) for the muon momentum dis- 
tribution. Correspondingly, we find r-appearance from the QE reaction as 
shown in Fig. |3](b). The event rate is 4.43 events for 414 days. 

The zenith angle distributions from the QE reaction are shown for fi (with 
and without oscillation) in Fig. |](a) and for r (from oscillation) in Fig. f|(b). 
As expected, all of the r produced are up-going. 

We do not discuss the details of DIS production of r because, as will be 
discussed later, such events are hard to identify. Nonetheless, we note that 
the event rate for DIS production is comparable to that of QE production. 
Therefore the current SuperKamiokande data set must contain about 10 r 
events for the oscillation parameters (|9|). 
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Figure 3: (a) fi and (b) r production rates from QE reactions at Super- 
Kamiokande for 414 days for lepton momentum greater than 1.33 GeV 
(Multi-GeV). The dashed histogram in (a) and the solid histogram in (b) 
are the case with neutrino oscillation, with parameters given in Eq. 
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Figure 4: The zenith angle distributions in (a) \x and (b) r production from 
QE reactions at Super-Kamiokande for 414 days in the multi-GeV region. 
The dashed histogram in (a) and the solid histogram in (b) are the case with 
neutrino oscillation, with parameters given in Eq. 
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3 Backgrounds 



The DIS production of r leads to multi-ring events. However, the character- 
istics of these events make it difficult to identify the r. The main difference 
from DIS CC production of e and fi is the absence of an energetic charged 
lepton, which is difficult to establish in multi-ring events, especially for e. 
The DIS NC events do not have charged leptons and are even more similar 
to DIS r events. Furthermore, the DIS cross sections for r production in 
the relevant range are smaller than those for the corresponding \x production 
by an order of magnitude. It appears to be difficult to separate the DIS r 
events from the DIS \x events, and we do not consider DIS r production in 
this letter. 

Instead, we focus on the QE production of r. Apart from the possible 
ring of the recoil nucleon, the entire event is caused by the decay products 
of the r. The purely leptonic decay modes, r — * u T uil for / = e, /i, have 
branching fractions of 17.8%, and 17.4%, respectively. These events would 
appear the same as the QE production of e or /i, and hence cannot be used 
as the signal. 

The single pion decay mode r ± — > has a branching fraction of 11.3%. 
The physics background to this signal is from the single pion production from 
neutral current (NC) interaction of the neutrinos with nucleons, 

Uin^> uipn , (10) 
viV^ Vimc, (11) 

where ^ can be any of the active neutrino species. In the energy range of 
our interest, the NC production of a single pion is dominated by the produc- 
tion of baryon resonances with their subsequent decays into a nucleon and a 
pion. Rein and Sehgal (RS) ]14[ attempted the calculation of the resonance 
contributions, assuming a simple harmonic oscillator model of the baryon 
resonances to calculate the transition amplitudes. Even though the model 
itself is somewhat questionable, their calculations show reasonable agreement 
with data within about 20% for charged current (CC) single pion production 



[15]. Therefore, we can rely on RS calculations on neutral current single 
pion production cross sections at 20% level. The cross sections for processes 
(HOJIl) are almost flat above E v > 3 GeV, both of them at 1.0 fb 0]. Un- 



fortunately, this is larger than the expected size of the signal cross sections 
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Table 1: The single rho production cross sections in fb using the resonance 





E u = 2 GeV E v = 20 GeV 


up — > vnp + 
vn — > vpp~ 


0.11 - 0.13 0.20 - 0.50 
0.074 - 0.14 0.26 - 0.57 



Eqs. (|7]||) once multiplied with the tau branching fraction to single pion of 
11.3%. This process also suffers from the "fake" background of charged cur- 
rent muon production where the muon is identified as a charged pion. This 
cross section is an order of magnitude larger, and is a continuum in muon 
energy with no features. Even though the SuperKamiokande Collaboration 
is working on the discrimination of charged pions from muons, which appears 
feasible based on their beam-based study at KEK [IB], identification of r in 
the pion mode appears quite difficult. 

The r — > v T p~ mode has a much larger branching fraction of 25.2%. We 
find this mode much more promising than the pion mode. Unfortunately, the 
total cross sections of the physics background, the NC single-rho production, 
are neither measured nor calculated. We therefore extend the RS calculations 
to rho production by summing over the resonance production cross sections 
given in their TABLE IV multiplied by the branching fractions into rho from 
Particle Data Group fl7fl . Our calculations neglect the interference between 
resonances as well as the continuum contribution, but we have verified that 
our naive method gives cross sections for the single pion production which 
agrees with their full calculations within about 10%. The larger uncertainty, 
however, results from the lack of data on the rho branching fractions of 
baryon resonances. We therefore give conservative and optimistic estimates 
of the NC single-rho production cross sections by taking the upper (up to 
100% if not measured) and lower (down to 0% if not measured) ends of 
the rho branching fractions, respectively. The results of our calculations 
are shown in the Table |l|. The uncertainties are unfortunately rather large. 
It is encouraging, however, that these cross sections are up to an order of 
magnitude lower than the signal cross sections of QE r production. 

There are other "fake" backgrounds to the single rho production from r 
appearance. The most important one is the charged current event with 
a single neutral pion production, i.e. v^n — > p~p7i° or v^p — > p + mr , where 
the muon is misidentified as a charged pion and the invariant mass of the 
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charged pion and the neutral pion falls into the range of the rho mass. The 
relevant cross sections are 



oiy^n^ pTpiP) ~ 3 fb, (12) 
vipfl -> p + nn°) ~ 2 fb. (13) 

Comparing with the signal cross sections, Eqs. (|7|J|), and r — > v T p branching 
fraction of 25.2%, we would like to achieve a rejection factor on p^n pro- 
duction by an order of magnitude (keeping S/N > 3). Even though we do 
not know the realistic number for it, it is reasonable to expect that the com- 
bination of the mass cut on 7r ± 7r° within the m p range together with 7r ± //i ± 
discrimination being discussed |nj, this level of rejection can be achieved. 



The next class of background events to discuss is the deep inelastic events 
with p production. The inclusive production rate of p ± from neutral current 
reactions is not measured. However, there is a reported rate of inclusive 
production rate of p° from charged current reactions, which give p multiplicity 
of about 0.1 |TJ| in the energy range of our interest. Assuming the comparable 



p production rate from neutral current reactions, together with the total 
neutral current cross section of about 7 fb at E v ~ 5 GeV, the inclusive p 
production cross section would be about 0.7 fb. Since these events are likely 
to be associated with production of a few other pions, the requirement of 
single p production is expected to suppress this background further. Overall, 
the background rate from deep inelastic p production is much smaller than 
the signal. 

Another possible source of background is the neutral current events of 
multi-pion production (DIS) where only one neutral pion is produced and 
only one charged pion is above the Cerenkov threshold. The estimate of this 
background depends sensitively to the detector capability and is beyond the 
scope of this paper. It should be noted, however, that there is no particular 
reason for the neutral and charged pions to form an invariant mass in the rho 
mass range for this process, and the continuum DIS process itself is relatively 
low at this energy range. 

As we have seen above, the estimates of the background are quite crude 
at best, primarily because of the lack of relevant data at this energy region. 

The efficiency in reconstructing rho is crucial in the study of r appearance. 
The closest study we could find in the literature is the search for proton decay 
p — > p + v by Kamiokande experiment | I9]| . The efficiency was unfortunately 



low: 15%. We can only hope that the Super-Kamiokande can achieve higher 
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efficiency thanks to much better granularity of the detector. It is encouraging 
to quote the high efficiency of 54% quoted by the 1MB collaboration in their 
search for p — > p ± v [ pOj] . 

Recall that we discussed only the QE contribution to r production. There 
are, however, additional processes of r production with additional particles 
where they are below the Cerenkov threshold. This adds more events to the 
single t production. On the other hand, the same enhancement occurs also 
for the background. It is not clear what the exact situation is. 

Let us emphasize that the background cross section can be measured 
by the near detector in the K2K experiment, the long-baseline experiment 
from KEK to Super-Kamiokande. The near detector consists of a 1 kt water 
Cerenkov detector followed by a tracking device, a calorimeter and a muon 
chamber. With the water Cerenkov detector they can study the "single- 
p" events. Then the uncertainty in the estimated background level will be 
significantly reduced. 

4 Statistical Significance of r Appearance 

We estimate the statistical significance for a r appearance signal at Super- 
Kamiokande. Suppose that the actual total event rate for signal and back- 
ground p events is 5* and B events per year, and that a fraction e of these are 
reconstructed. The background events will be up-down symmetric, whereas 
the signal events are dominantly up-going. We assume a value for Am 2 such 
that all signal events are up-going. With N years of data there will be NeB/2 
down-going events and Ne(B/2 + S) up-going events. 

For quasi-elastic r production we have calculated S = 1.0 for sin 2 29 = 1; 
while for single p production we estimate B = 0.4, using the extension of the 
Rein-Sehgal resonance calculation discussed in the previous section. Hence 
with e = 0.5, and N = 4(8) we expect 0.4 (0.8) downgoing events and 
2.4 (4.8) up-going events. We would like to determine the probability that 
the up-going events are caused by a statistical fluctuation of the background; 
however, it is difficult to do this precisely because there are so few down-going 
events that the background is very poorly determined. Hence we will assume 
that the background has been determined from measurements in a neutrino 
beamline, such as at the near detector of the K2K experiment. In this case, 
for the observed excess of up-going events compared to down-going events to 
be due to a statistical fluctuation of the background requires a fluctuation 
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of y/2Ne S/y/B standard deviations. With S = 1 and B = 0.4, this gives a 
4.6<t effect for iVe = 4. As discussed earlier, there are other contributions to 
both S and B, so that it is probable that the S/y/B value used here, 1.6, is 
not correct. For S/y/B = (1.0, 1.6,2.5), the signal is a (2.8, 4.6, 7.1)a effect, 
for iVe = 4. Hence a significant result is only expected from a multi-year 
search with high p reconstruction efficiency The event rates are low, so that 
in practice the signal could be worse or better depending on whether the 
number of observed events happens to fluctuate above or below the mean 
expected number. 

5 Conclusion 

We have studied the prospect of detecting r appearance in Super-Kamiokande 
detector, from the oscillation of atmospheric to u T , and have identified 
r — > v T p^ as a promising detection mode. The event rate, however, is low: 
about 1.0 event per year from quasi-elastic production. There are possibly 
some additional contributions from non-QE processes where the other prod- 
ucts (pions) are below the Cerenkov threshold. We estimated rates for back- 
ground processes, and, even though the estimates have large uncertainties, 
none of the background processes appear larger than the signal. The S/N 
can range anywhere between ~ 1/3 to 1. This uncertainty will be reduced by 
direct measurements of the background at the K2K experiment. Although 
challenging, with 5 to 10 years of data, the detection of r appearance from 
atmospheric neutrinos may be possible, for example, with an observation of 
5 upward going events and 1 downward. 
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